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Abstract 

We present a quantitative semiclassical treatment of the effects of bi- 
furcations on the spectral rigidity and the spectral form factor of a Hamil- 
tonian quantum system defined by two coupled quartic oscillators, which 
on the classical level exhibits mixed phase space dynamics. We show that 
the signature of a pitchfork bifurcation is two-fold: Beside the known 
effect of an enhanced periodic orbit contribution due to its peculiar Ti- 
dependence at the bifurcation, we demonstrate that the orbit pair born at 
the bifurcation gives rise to distinct deviations from universality slightly 
above the bifurcation. This requires a semiclassical treatment beyond the 
so-called diagonal approximation. Our semiclassical predictions for both 
the coarse-grained density of states and the spectral rigidity, are in excel- 
lent agreement with corresponding quantum-mechanical results. 
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1 Introduction 



A prominent approach to the quest of "quantum chaos" involves spectral statis- 
tics to characterize the energy-level fluctuations in quantum systems and their 
interpretation in terms of the dynamics of the corresponding classical system. 
Classically integrable systems possess uncorrelated energy levels, described by 
a Poisson distribution [T], while the levels of classically chaotic quantum sys- 
tems exhibit strong local repulsion. This behaviour is conjectured to be the 
same as for the eigenvalues of ensembles of random matrices preserving certain 
general symmetries [2]. Spectral statistics has been investigated, for both inte- 
grable [21 |H[S] an d chaotic [HI El IE] systems, employing semiclassical (periodic 
orbit) approaches, which provide the closest link between classical and quan- 
tum properties. For the purely chaotic case, starting with Ref. [5], considerable 
progress has been recently made in understanding energy level correlations semi- 
classically beyond the so-called diagonal approximation [7] by means of classical 
correlations between (off-diagonal pairs) of periodic orbits [TU] . 

However, intcgrability and full chaoticity represent extreme situations which 
occur rather exceptionally. The most realistic physical situation is that of a sys- 
tem which is neither completely chaotic nor integrable, but whose phase space 
contains a mixture of stable orbits surrounded by regular islands and chaotic re- 
gions. One main feature and structuring clement of classical mixed phase space 
dynamics is the occurrence of bifurcations of periodic orbits upon variations 
of the energy or other parameters of the Hamiltonian. Moreover bifurcations 
lead to noticeable effects in the spectral statistics, because in semiclassical trace 
formulae for the density of states [Til Q~2] , contributions from periodic orbits at 
a bifurcation exhibit an enhanced weight, compared to that of isolated orbits. 
This has been demonstrated for the generalized cat map in Ref. [13] , where the 
semiclassical signature of a tangent bifurcation was studied on the level of the 
diagonal approximation. 

More generally, in Ref. [T4J, [15] a semiclassical approach was developed for 
the moments of the level counting function in the presence of several compet- 
ing generic bifurcations. It was suggested that these moments diverge with a 
universal "twinkling exponent" in the semiclassical limit fi — > 0. 

In the present paper we investigate the role of pitchfork bifurcations on the 
spectral statistics in Hamiltonian systems that are closer to a realistic physical 
situation than the maps considered so far. We show that bifurcations of short 
orbits has a considerable effect on the spectral rigidity and the spectral form 
factor, respectively, even in the almost chaotic case. As a standard system with 
mixed classical dynamics, we choose the Hamiltonian of two coupled quartic 
oscillators. Its relevant classical bifurcation characteristics is summarized in 
Sec. [2] In Sec. [3] we present a detailed semiclassical analysis including a com- 
parison with quantum results for the (smoothed) density of states for different 
symmetry classes, as a prerequisite for the treatment of spectral correlations 
in Sec. [U There we quantitatively analyze deviations of the spectral rigidity 
from universality employing uniform approximations to derive the semiclassical 
periodic orbit weights at the bifurcation. We show, in particular, that pairs of 
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orbits (with an action difference smaller than Planck's constant h), born at a 
pitchfork bifurcation, yield important non-diagonal contributions to the spectral 
form factor and rigidity. The deviations from the quantum chaotic universality 
are found to be most significant after, rather than at the bifurcation. 

2 The quartic oscillator Hamiltonian 

As a representative system we investigate the coupled quartic oscillator (QO) 
in two dimensions. Its Hamiltonian reads: 

H{x,y, Px , Pv ) = -(p 2 x +P 2 y ) + - (x 4 + y 4 ) + ! x 2 y 2 . (1) 

It has been extensively studied both classically, semiclassically and quantum- 
mechanically [HI [T71 [TSl [THl HO] , as a smooth potential model exhibiting the 
transition from integrability to chaotic behaviour. Here we summarize the main 
classical features relevant for the subsequent semiclassical treatment. Since the 
Hamiltonian fl} is homogeneous, its classical dynamics at different energies E 
can be related to each other by a simple scaling of coordinates, momenta and 
time. All actions scale with energy E as E 3 ^ 4 so that the semiclassical limit can 
be unambiguously taken as E — > oo. 

After scaling out the energy the parameter a in Eq. fl]) solely determines the 
dynamics. The system is integrable for a — 0, 1, and 3. For a = 9, it is almost 
completely chaotic: we could not locate any stable periodic orbit with a period 
up to about four times that of the shortest orbits. For values a > 9 the regular 
fraction of the phase space keeps oscillating with a decreasing amplitude. 

The QO in Eq. ([1]) possesses periodic straight-line librational orbits along 
both axes which we label by A. They undergo stability oscillations under the 
variation of a. Infinite cascades of new periodic orbits bifurcate from the A 
orbits and their repetitions. The motion of the A libration can be given analyt- 
ically in terms of Lame functions [TH1 HH] • The trace of its stability matrix M 
(see [TTJ [TJ] for its definition) as a function of a is known analytically [3T] : 

TrM(a) = 4cos (|\/T+8a) +2. (2) 

Isochronous pitchfork bifurcations of the A orbit (which are non-generic due to 
the discrete symmetries of the system) take place when TrM=+2, i.e., for 

a = a n = — n (n + 1) , n = 0, 3, 4, 5, . . . . (3) 

(For a.\ = 1 and a2 = 3, where the system is integrable, the A orbit is member 
of a degenerate family and does not bifurcate. See also [HI QI5] for more details 
about the periodic orbits of this system.) 

In Fig. [T] we show TrM(a) for the primitive A orbit and the new orbits 
born at its bifurcations at a n with n — 3 to 7. These orbits are alternatingly 
stable or unstable rotational (R CT ) and librational orbits (L^) with a classical 
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degeneracy of 2 due to the symmetries (cf. Ref. [E].) In our numerical case 
studies below, we shall focus on the bifurcation at a = = 10 where the orbit 
L6 is born. Note that at each second bifurcation (n = 3, 5, . . . ) a new stable 
orbit (R5, L7, . . . ) is born, so that stable orbits exist on either side of these 
bifurcations. At the other bifurcations (n = 4, 6, . . .), on the other hand, the 
new orbits (L 6 , R 8 , . . . ) are unstable, and just before these bifurcations, the A 
orbit is also unstable. This explains the oscillating regularity of the phase space 
and the fact that, even in the limit a — > 00, there always exist regions with 
stable orbits [2D]). 




Figure 1: Trace of the stability matrix M as a function of a, Eq. (fl]), for the 
primitive A orbit (solid line) and the new orbits born at its bifurcations (dashed 
lines) at a = 6, 10, 15, 21, and 28. Subscripts denote the Maslov indices Oj (see 
Sec.©. 

The potential in Eq. ([T]) is invariant under the symmetry operations that 
conform the point group symmetry C4y, which has four one-dimensional ir- 
reducible representations and one (doubly-degenerate) two-dimensional repre- 
sentation. Due to the C<±v symmetry, the full eigenvalue spectrum would not 
exhibit any universal statistics. For an appropriate study of the spectral statis- 
tics, each symmetry class must therefore be treated separately. We shall study 
mainly the representation corresponding to cigenfunctions which are symmetric 
under the operations x — > —x, y — > — y and x — > y, which we call EES. This 
representation is easier to handle semiclassically, because all its characters are 
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equal to unity. 

For the numerical calculation of the quantum mechanical eigenenergies we 
follow the procedure outlined in Ref. [55J . We diagonalise the Hamiltonian using 
a basis of symmetry-adapted linear combinations of harmonic oscillator states: 

\n x ,n v ) m = -^=(|ni,n 2 ) ± |n 2 ,ni)) i ( 4 ) 

where the sign and the parity of n\ and ri2 depends on the representation. Since 
the independent symmetry-reduced blocks of the Hamiltonian matrix in this 
basis are banded, we can obtain up to tenthousand well-converged eigenvalues, 
allowing for significant statistics. 



3 Semiclassical density of states for discrete sym- 
metries 

Periodic orbit theory yields the semiclassical spectral density as 

g(E) = d(E) + Sg(E) (5) 

where the smooth part d{E) is given by the (extended) Thomas-Fermi model 
(cf. Chap. 4 in Ref. [12]), and the oscillating contribution is given by a trace 
formula which, to leading order in \/Ti, has the following form: 



W) = ^iE ME) cos 



Sj(E) 7T 

IT - 2°* 



(6) 



The sum is over all periodic orbits j (which form families with degenerate actions 
in the presence of continuous symmetries). Sj(E) — §. pdq is the action integral 
along a periodic orbit and (Tj a geometrical phase factor (usually called Maslov 
index). The amplitudes Aj(E) and the power of U in Eq. ([6]) depend on the 
presence of continuous symmetries. For systems without continuous symmetries, 
where all orbits are isolated in phase space, one has \i = 0, and the amplitudes 
Aj(E) were given by Gutzwiller [53] in terms of their stability matrices Mj(E) 
and periods Tj(E) — dSj(E) /dE. When an isolated periodic orbit undergoes 
a bifurcation at an energy Eq, its amplitude in the Gutzwiller trace formula 
diverges and uniform approximations must be developed [24] to obtain a finite 
Aj(Eo); in this case one finds < fi < 1/2, the precise value of n depending 
on the generic type of the bifurcation (cf. also Ref. [I25J. For fully integrable 
systems, \i = //2, where / is the degree of degeneracy of the most degenerate 
orbit families; the amplitudes were derived by Strutinsky and Magner |26j for 
specific cases and by Berry and Tabor [27] for general integrable systems (cf. 
also Sec. 13.11 below). For non-integrable systems with continuous symmetries, 
further results were obtained by Creagh and Littlejohn [55], who also derived a 
Berry- Tabor- like trace fromula for integrable systems. 
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In the presence of discrete symmetries it is necessary to define partial den- 
sities of states corresponding to the subspectra of each irreducible represen- 
tation of the symmetry group. For systems with isolated orbits, the corre- 
sponding symmetry-reduced semiclassical trace formulae have been derived in 
Refs. [2S] GUI EH]; we shall discuss and use them in Sec. 13.21 

For practical purposes, it is useful to coarse-grain the density of states by con- 
volution with a normalized Gaussian exp[— (E/'y) 2 ]/(y/n r y). Hence, we replace 
the quantum density of states d(E) = J2 n $(E — E n ) by the "coarse-grained" 
density of states 



exp 



(E - E n f 



r 



(7) 



whereby the smoothing width 7 defines the energy resolution at which one 
wishes to study the spectrum. The correspondingly averaged semiclassical level 
density becomes, to leading order in h (see, e.g., Ref. [12J, 



5 9l {E) 



2K 



(8) 



Hence, long orbits are exponentially suppressed which avoids convergence prob- 
lems for not too small values of 7. 



3.1 Integrable Systems 

For integrable systems with / degrees of freedom, it is useful to work with 
action-angle variables (I, (pi), with each set of actions I = {Ii, ...,//} defining 
a phase-space torus [32]. The Hamiltonian can be transformed to H(T) = E, 
and the frequencies d(fi/dt = u = {u>i, . . .,oj/} on the torus I are given by 
<jj(T) = V-ff(I). Assuming smooth boundaries, the Einstein-Brillouin-Keller 
(EBK) quantization [33 



T j (n j ) = h(n j + 1/2) 



nj = 0,1,2,,. 



defines a set of / quantum numbers n = (tij., 
into E = H(T), the EBK spectrum reads 



,71/). 



i = 1, • • • , / , (9) 

Upon inserting Eq. ^ 



tpEBK rpEBK 
-^n = i 'n 1) ...,n / 



H(h( ni ),...,I f (n f )) 



(10) 



Berry and Tabor [57] started from the density of states in terms of the E^ BK 
and converted it, by means of Poisson summation, into a semiclassical trace 
formula of the type of Eq. ©. 

The EBK quantization of the integrable QO, Eq. ([T]), with a = has been 
performed in Ref. [19j : we quote here those results which are relevant for our 
present application. The EBK spectrum becomes 



E 



EBK = \( 



{n x ,n y =0,1,2,...) 



(11) 



G 



where K = K(k) is the complete elliptic integral of first kind with modulus 
K = l/v2' Since the Hamiltonian ([1]) is separable for a — 0, we can write 
E„ B n — En BK + E EBK ■ The separate one-dimensional densities of states, 



9j (E) =J2 6 { E - K BK ) . 0' = *> y) (12) 

rij— 

which are identical due to the symmetry, become after Poisson summation 

9]{E) = £ ( ~ 1)fej cos M^)A] > o* = *,v) (13) 
fc 3 =i 

corresponding to the Gutzwillcr trace formula for a one-dimensional system. 
Here 

S A {E) = Ut{AEf/\ (14) 

is the action of the primitive A orbit and Ta(E) — &Sa{E)/&E its period. The 
total density of states of the full two-dimensional system can then be written as 
a convolution integral of the one-dimensional densities: 

g(E)= / g x (E-E')g y (E')dE'. (15) 
Jo 

The asymptotic evaluation [33] of this integral in the limit fi — > yields for the 
oscillating part 



Sg(E) - 2(^) 2 (4E)iy y (-l)^+ fc « ^ 5 cos 



(16) 



The double sum in the first line above contains the contributions from the 
standard stationary-phase evaluation of the integral. It corresponds exactly to 
the Berry- Tabor trace formula [27j , whereby the two numbers k x , k y label the 
rational tori corresponding to the simply degenerate families of periodic orbits 
with two-dimensional motion. The actions of these rational tori are given by 

S kxkv {E) = S A {E){kl + kl) 1 '\ (17) 

The term in the second line of Eq. (|16p arises from the boundaries of the 
integral (|15p , corresponding to the A orbits which are one-dimensional librations 
with all energy in either x (E' — 0) or y direction (E' = E). Note that the 
amplitude of this term involves a prefactor Ti~ 5 ^ 4 . This is due to the fact that 
the A orbit undergoes a pitchfork bifurcation at a = corresponding to n = 
in Eq. ((3j. (The orbits L3 born at this bifurcation exist only for a < 0.) 
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In Ref. [12], identically the same result (flB]) was obtained, whereby the local 
uniform approximation [24] for the contribution of the bifurcating A orbit was 
employed [55] . 

In the upper panel of Fig. [5] we compare the semiclassical density of states, 
Eq. (|16|) (dashed line), with the corresponding quantum- mechanical one (solid 
line), both coarse-grained with a Gaussian average with width 7=1. We find 
perfect agreement up to very high energies. 
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Figure 2: Upper panel: total density of states for a = coarse-grained by a 
Gaussian with width 7 = 1. Lower panel: symmetry-reduced density of states 
for the representation EES, see text. Solid line: quantum result, dashed line: 
semiclassical result, Eq. (|16l) . 

We now calculate the symmetry-reduced densities of states by restricting our- 
selves to the subspectra, E BBK , of a given irreducible representation. Hereby 
we can relate the parities of the quantum numbers to the symmetries of the 
irreducible representations. Thus, we restrict n x and n y to be even or odd, ac- 
cording to a given representation. For example, let us take the one-dimensional 
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irreducible representation EES. This corresponds to taking n y < n x with n 



even. Then the partial density of states can be calculated as a convolution 

r-E 



Sg^(E) 



of the one-dimensional densities gf(E) defined as in Eq. (JT3J) , except that only 
the terms with even rij are included in the sum. The asymptotic evaluation of 
the convolution integral leads to 



Sg^(E) 



K 

itTi 



1 (K)t 
2i (xh)i 



1 



2h 



Skxk y {E)--(k x + k y )- 



(AE)-T-ey — 



2h y ' 2 8 



(18) 



Again, the first term above corresponds to the Berry- Tabor result for the ratio- 
nal tori, and the second term comes from the bifurcating A orbit. 

In the lower panel of Fig. [5] we compare the semiclassical and quantum- 
mechanical density of states, 5g EES (E), coarse-grained with a Gaussian average 
with width 7 = 1. Again the agreement is nearly perfect. 



3.2 Isolated orbits 

The symmetry-reduced densities of states for isolated orbits have been derived in 
Ref. [29l [30] by projecting the semiclassical Green function onto the irreducible 
representations and reducing the classical dynamics to the fundamental domain 
which is the smallest part of the phase space which tesselates the whole space 
under application of the allowed symmetry operations. After this procedure one 
obtains the reduced density of states in the irreducible representation m 

Here d m is the dimension and Xm(g) the character of the symmetry operator g 
in the irreducible representation m. The bars in Eq. (TT9|) indicate that actions, 
periods, stability matrices and Maslov indices are calculated in the fundamental 
domain, while g\ is the operator that relates the r-th repetition of the reduced 
orbit I with its original lifted into the the whole phase space. \Ki \ is the order of 
the group Ki which leaves every point of the orbit I invariant. By the definition 
of the fundamental domain, this is the identity for orbits that stay in the interior 
of the fundamental domain, while there can exist more than one operation for 
orbits that lie on the boundaries. The matrix D; is block-diagonal in coordinates 
with blocks given by d(gq)/dq with g £ K. This matrix is again the identity for 
interior orbits, but can be different for boundary orbits. 

It is usually easier to solve the equations of motion in the whole space than in 
the fundamental domain, where one has hard- wall reflections. Given the classical 
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quantities for the total space, the task is then to find their reduced counterparts 
(marked with bars in (|19p ). Take a Hamiltonian of the form H(p, r) = p 2 /2m + 
V(r) which is invariant under the point-group symmetry G. Suppose that the 
subgroup H leaves the / orbit invariant (not pointwise), then the I orbit can be 
divided into \H\ copies related by symmetry [35]. There will be |G|/|7J| copies 
of the orbit in the full phase space (if we consider time reversal, then there are 
2|G|/|JJ| copies of orbits without time-reversal symmetry). Therefore the lifted 
orbit should be equivalent to the \H\ = r-th repetition of the reduced orbit (or 
to the | H | /2-th repetition for time-asymmetric orbits, which become librating 
orbits in the fundamental domain, and the |ff|/|if|-th repetition for boundary 
orbits). Hence, all the classical quantities should be inter-related as 



Si(E) =rSi{E), T l (E)=rT l {E), ai=fa u Mj = M[ , (20) 



since they are invariant under point transformations. The only difficulty remains 
to find out which of the roots of M; must be taken. E.g., for = 2 we have 

Mj = M; . Thus, if the eigenvalues of M; are e ±Ul , those of M; can be ±e ±Ul/2 . 
On the other hand, we know that for two-dimensional Hamiltonian systems, 
hyperbolic orbits always have even Maslov indices, while elliptic and inverse- 
hyperbolic orbits always have odd Maslov indices [37]. We have observed that 
this rule can be reversed in the fundamental domain. 

This is illustrated in Fig. [3J for the case of a single reflection symmetry with 
respect to the x axis. Then the fundamental domain is the upper plane (y > 0) . 
We have calculated the Maslov index W using the method of Creagh et al. [38] 
(as explained in Ref. [12] . App. D) and verified that it is, indeed, either the 
same as for the lifted orbit for orbits without this symmetry, or half of it for 
orbits with reflection symmetry. However, the sign of the eigenvalues did not 
follow Sugita's rule [37]. This rule can, however, be applied to W — mod(R,2), 
where R indicates the number of hard-wall reflections at the boundaries of the 
fundamental domain. Thus, if this number is odd, the rule is reversed. 

We have calculated the reduced density of states (fT9"]) for the representation 
EES in the QO at a — 9. The result is shown in Fig. @] for Gaussian smoothing 
with width 7 = 4. A considerable agreement between the semiclassical (dotted 
line) and the quantum-mechanical result (solid line) is achieved. 



To study the effect of pitchfork bifurcations on the spectral statistics we consider 
the spectral rigidity or stiffness, A [35]. It is defined as the local average of the 
mean-square deviation of the staircase function N(E) from its best-fit straight 
line over an energy range corresponding to L states with mean level spacing d: 



4 Spectral Rigidity 




de [N(E + e) - A 




(21) 
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Figure 3: Calculation of Maslov indices for some reduced orbits of the QO at 
a = 9, considering only the reflexion symmetry at the x axis. Left panels: 
reduced orbits in (x,y) plane. Right panels: evaluation of the Maslov index a 
which corresponds to the winding number of the complex number C (t) over one 
period (cf. [38). Top panels: librational orbit B2 along the diagonal. Here the 
length of the reduced orbit is the same as that of the lifted orbit, and their 
Maslov indices are equal. Center panels: orbit R 4 . Here the reduced orbit is 
half of the lifted orbit and its Maslov index is ef = 2 (i.e., half of the total a) but 
TrM is negative in spite of the even Maslov index. Bottom panels: orbit A§. 
The reduced orbit is again half of the total orbit, and so is the Maslov index. 
But Tr M is positive in spite of the odd Maslov index. 
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Figure 4: Reduced density of states for the representation EES in the QO at 
a — 9 after Gaussian averaging with width 7 = 4. The solid line shows the 
quantum result and the dotted line the semiclassical result using Eq. (TT9"|) . 

The quantity A(L) measures spectral correlations over energy distances of order 
L. For an uncorrelated Poisson spectrum the universal prediction is 

A Polsson (£) =L/15, (22) 

while for a chaotic system it is approximately given by 

A RMT (L) = JLlogL-D, (23) 

where D is a constant, (5 = 1 for systems without time reversal symmetry (GUE 
statistics) and (3 = 2 for systems with time reversal symmetry (GOE statis- 
tics). This universal behaviour has been observed up to correlation lengths 
L < i max = 2irhd/T- m i n , where T m i n is the period of the shortest orbit. In 
Fig. [5] we show the numerical results for the quartic oscillator in the integrable 
and almost chaotic regime, compared with the corresponding predictions, Eqs. 
l[2"2l |2"3"1) . The L range, in which the numerical data coincide with the universal 
predictions, increases with increasing energy, i.e., by approaching the semiclas- 
sical limit. 

For a mixed system it was conjectured that the statistics will be a super- 
position of Poisson and random matrix contributions [401 141] ■ parameterized 
as 

A(L) « A Poisson ((l - q)L) + A RMT (<ji), (24) 

were q is the irregularity fraction of the system (i.e., the fraction of the phase 
space corresponding to the chaotic sea). Since both statistics are monotonously 
increasing functions, we expect that the more regular the system is, the larger 
is the rigidity. 
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Figure 5: Rigidity for a = (integrable case) and a — 9 (almost chaotic 
case). With increasing (unfolded) energy E the numerical data converge to 
the universal Poisson (left panel) and random matrix predictions (right panel) 
marked as full lines. 



4.1 Semiclassical theory for the rigidity 

The semiclassical theory for the rigidity was developed in Ref. [7], for the two 
limiting cases of complete chaoticity and full regularity (integrability) . The 
procedure is the following: By energy integration of the density of states, Eq. ©, 
one obtains an expression for the number of states. By inserting this expression 
into the definition of the rigidity one finds 



A(L) 



1 



EE' 



2ft 2 " \^ V T i T * 

7 k J 



where Tj — dSj/dE are the periods, 



LT) 

2hd 



(Sj - S k ) + -fa - crfe) 



T 4 



G{y 3: y k )), (25) 



(26) 



and 



G(x,y) 
F(x) 



F(x -y)- F(x)F(y) - 3F'(x)F'(y), 
- sin x = j (x) . 



(27) 
(28) 



The main contributions come from pairs of orbits whose action difference is 
smaller than h, so that yj can be chosen to be equal to y k in the argument of 
G: 



A(L) = 



1 



AjA k 



j k 



exp 



Sk 



7r ' i 



Ok) 



2ft 2/i 

where y~jk = \{yj + yk) and g{x) — G(x,x). The function g(y) (see Fig. [6J 
selects the orbits that contribute to the double sum. If L « L max then g(y) 
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is almost unity only for long orbits, while for L > L max the function is almost 
unity for all y, and the most important contributions to A(L) come from short 
orbits due to the factor 1/T 2 . Since we are interested in studying the effects of 
a bifurcation of one of the shortest orbits, we are going to concentrate on the 
saturation behaviour, which corresponds basically to the first moment of the 
staircase function. 

The rigidity can be written in terms of the spectral form factor K(t) (the 
Fourier transform of the autocorrelation fuction) as 



1 f K(t) 
A W = ^ J -y^Lr)dr, 



(30) 



with r = T/2-Khd and K(r) = A JZo( d ( E + u/2)d{E - lu/2)) e^^da;) . 

\ / At 

A local time average At has to be performed in order to obtain a self-averaging 
form factor. 

The corresponding semiclassical expression for the form factor, analogous to 
Eq. ([291), is 



K{t,E) 



T 2 

1 H 



2^ 



Cfe) 



AE 



(31) 

Tfc). The width of the delta- function is due to the time 



where T jk = |(Tj 
average At. 

As expressed in Eqs. ([251 l3~Tj) . the rigidity and the spectral form factor are 
determined by a double sum over pairs of periodic orbits. The semiclassical 
limit h — > means that the typical classical actions of these paths are very 
large compared with h, so that the energy average will strongly suppress the 
contributions of most pairs of orbits. The first approximation is to consider that 
only orbits paired with themselves (j — i) or with their time-reserved partners 
(j = *) §i ve a contribution, which is known as the "diagonal approximation" 
0. 




Figure 6: Window function g(y), see text. 
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For the QO at a = 0, the tori amplitudes A^ x ,ky are given by 



A 



kx,ky 



3/2 



(k x + ky)^^ 



(32) 



for the irreducible representation EES. For integrable systems the contribution 
of the non-diagonal terms j : k in the sum (|25p will vanish after averaging 
owing to destructive interference. For this system, due to the degeneracy in 
the actions, the orbits that contribute to the double sum are those that satisfy 
n x + n y = n 'x + n 'y- Inserting the amplitudes for the tori and summing only 
over terms with the same actions we have 



A(L) 



(415) 



3/4 



2 4 tt 3 ?jK 



E 



h h 

■g(yk x ,k y ) 



n x n y o h - U 



(33) 



where Ik = + ky) 1 / 4 , l n — (n 4 + n 4 ) 1 / 4 , and 6 is the Kronecker delta. 

With this expression, we can reproduce very well the statistics semiclassi- 
cally, as is shown in Fig. [71 




e 1 1 1 1 

50 100 150 200 

L 

Figure 7: Semiclassical (dashed lines) and quantum results (solid lines) for the 
spectral rigidity A(L) for different values of E. 

We have neglected here the contribution of the A orbit which undergoes 
a pitchfork bifurcation. We have checked that its contribution is negligible, 
since its amplitude in the PO expansion goes like ?i~ 5 / 4 (a power one quarter 
larger than an isolated orbit) compared with that (h~ 3 ^ 2 ) of the tori. For the 
saturation we can take g(yk x ,k y ) — 1- Then the energy dependence of the 
saturation value Aqo goes like E 3 / 4 , as seen from Eq. (|33[) . 
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In the left panel of Fig. [5] we depict the saturation value obtained from 
the quantum spectrum (dots), which is well reproduced by the semiclassical 
prediction (solid line). For the chaotic case, RMT gives a saturation value A^ 



< 2 
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Figure 8: Saturation value A^, plotted versus unfolded energy E. The dots 
mark the quantum results. Left panel: integrable case (a — 0); the solid line 
shows the semiclassical prediction. Right panel: Almost chaotic case (a = 9). 
Here the solid line represents the GOE prediction. 



that behaves as log(l/?i) which is obtained by replacing the form factor by its 
GOE prediction in Eq. ([50]) . Though the exact saturation value is not universal, 
since it depends on the lower integration limit r m j n , its Ti dependence is. In the 
right panel of Fig. [5] we compare the quantum result with the GOE prediction 
evaluated for our value of T m j n for a = 9. 



4.2 Bifurcation effects in the rigidity 

It has been discussed in Ref. [13] that additional contributions to the long-range 
spectral correlations may arise from bifurcations of periodic orbits, and that this 
effect can be reproduced semiclassically. The authors of Ref. [T3] investigated 
the cat map at a tangent bifurcation, and found that the number variance of 
the counting function shows a "lift off" reaching a much higher value than in 
the normal chaotic situation. We report here similar findings for the rigidity of 
the QO Hamiltonian for values of a near the pitchfork bifucations of the A orbit 
at a n . Moreover, we find that the increase of the saturation value Aoo becomes 
even larger slightly above the bifurcations. This is illustrated in Fig. [9] In the 
left panel we show the rigidity A(L) for four values of a around a — = 10 
where such a bifurcation occurs. The rigidity at 04 = 10 exhibits a slightly 
larger saturation than at a — 9 ("lift off"). However, the increase is even much 
more noticeable at a = 10.5. Then the saturation goes down again for a = 11, 
even though the system is more regular than at a — 10.5 [32] . 

The energy dependence of A m is shown in the right panel of Fig. [5J We 
see that this effect exists over a large region of energies. As depicted in Fig. [TT] 
the phase space looks completely chaotic at the bifurcation at a = 10; without 
knowledge of the bifurcation one would expect an almost universal behaviour. 
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Figure 9: Left panel: spectral rigidity for a = 9 (crosses), a = = 10 (filled 
squares), a = 10.5 (circles), and a = 11 (triangles) for E = 4000. Right panel: 
saturation value Aoo versus E before and after the bifurcation at = 10. 
Although the phase space is barely affected, the saturation at a = 10.5 is much 
larger than the saturation at a = 11. 




Figure 10: Left panel: Form factor at a = 9, 10 and 10.5 compared with RMT 
(line). Note the strong peak at r = ta — 0.015 (indicated by an arrow) coming 
from the bifurcating orbit. Right panel: Form factor at a — 9 (crosses), a = 10 
(squares), a = 10.5 (circles), and a = 11 (triangles) in a zoomed region around 
t ~ ta- For a — 10.5: the amplitude of the peak is clearly larger than at the 
bifurcation. 
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a = 10.5 



Figure 11: Poincare surfaces of section for the QO near 04. At a — 9 and at 
the bifurcation point a = 10 the phase space looks competely chaotic. A new 
stable island appears at a — 10.5, which is slightly larger at a = 11. 



Above the bifurcation, a tiny regular island is seen at the center, which arises 
from orbit A7 that became stable. The island is slightly larger at a = 11 than 
at a = 10.5 (see Fig. ITT]) . 

Equivalently, in Fig. [TU] we show the effect in the spectral form factor. In 
the left panel we show K{t) at a = 9,10 and 11. The results are consistent 
with the GOE prediction for almost all times, but we see a very large peak at a 
time that corresponds to the period of the libration orbit, ta- This is consistent 
with the results of 13J. However, the enhancement is even more noticeable at 
a = 10.5 (right panel). 

The exact calculation of the semiclassical rigidity for the QO in the chaotic 
regime is numerically impossible, since this would require an infinite number 
of periodic orbits, and there is no analytical way to calculate them. To repro- 
duce the quantum result semiclassically, we calculate the coarse-grained reduced 
density of states, defined analogously to Eq. ([5]) by 



WE) 



1 



\Ki\' 



E 



Xrn(9l) 
\Wt-Th\ 



Si(E) 



(34) 



The longer orbits will be exponentially suppressed assuring convergence, but, 
at the same time, affecting the universality. However, for the study of the 
saturation properties of A(L) as a probe for bifurcation effects, the information 
of the shorter orbits should be sufficient. 

Consistently we also coarse-grain the quantum stair-case function, defining 



N^E) 



IE 



1-erf 



7 



(35) 



Inserting N y (E) into Eq. ([2Tj) , we obtain a "smoothed" rigidity A 7 of the coarse- 
grained density of states. We find that even for relatively large values of 7, the 
bifurcation effects described above are still clearly visible, as shown in Fig. [T2J 

We are now equipped to calculate the saturation property of the smooth 
rigidity A 7 (L) semiclassically, taking into account the bifurcation of the A orbit 
at «4 = 10. Its contribution to the total density of states, together with that of 
the L6 orbits born at the bifurcation, to the total density of states is given in 
the "global" uniform approximation of Ref. (with a = +1, a < 0, <ti = — 1 
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and v = <JA & ,r = 6r for the present case). It reads 



SgT +L (E) = $le\ 
irn 



x |A 



ttAS" 



2h 



1/2 



exp I — S — i3rir — i— I x 
a 4 



^ /4 (^) e -, +J _ 1/4 (l^ 



(36) 



Here AA = A L /2 - A A /V2, A = A L /2 + A A /V^, S={S L + S A )/2 and AS* = 
(Sl — Sa)/2, where Aj(E) and Sj(E) are the Gutzwiller amplitudes and actions 
of the isolated A and L orbits, respectively, away from the bifurcation, r is their 
repetition number, and oi = sign(a — a^). At the bifurcation (a — = 10), 
the local uniform approximation becomes 



Sg. 



2TrV2^h 5/4 \a\ 1 / 4 r 3 /' 



■ cos 



Sa 



3r7r 



(37) 



Here Ta{E) is the period of the primitive A orbit, and a is a normal form 
parameter which we determined numerically from the local expansion given in 
Eq. (gl]) below (cf. also Ref. pH [25]). 

In this way we can reproduce the quantum mechanical results near the bi- 
furcation semiclassically, as demonstrated in Fig. [13J Further analysis showed 
that amplitudes and actions of most of the orbits do barely change, and the 
higher saturation for the smooth rigidity was mainly caused by the bifurcation. 

Considering the rigidity without smoothing, we now assume that the contri- 
bution of the long orbits corresponds to and can be replaced by the universal 
RMT prediction, so that the differences in the saturations arise basically from 
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Figure 13: Smoothed rigidity for a — 9, a = = 10 and a = 10.5, ob- 
tained for 7 = 2. The solid (dashed) curves represent the quantum mechanical 
(semiclassical) results. 



the A and L orbits. Hence, we approximate the saturation value of A by 

Aco(£) 



pGOE , ^A,L 



A 



GOE 



2 \ ^ T,T k 

j,k=A,L 3 K 



■ cos 



(38) 



At the bifurcation, the second term corresponds to the diagonal contribution of 
(|3"7|) , so that 

r 2 (i/4) 



, A,L 



87r 3 |a|V2nV2 ' 



and behaves like 



ioo oc log(l/?i) + 



(39) 



(40) 



In the neighborhood of the bifurcation, i.e., when the action difference IAS*! is 
smaller than h, we can expand the actions and amplitudes around a = (cf. 
Ref. [55]): 



AS" 



Sa — S L 



A. 



0(6 3 ), 

[l + 0(e)] ; 



2 4a 

l± t Al = Ia 

where e = c (a — 014). Up to first order in e this yields 

T4 



(41) 
(42) 
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(43) 
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Inserting this into the saturation value of the rigidity we obtain 



A A ' L « + e _J_ + c * r2 ( 3 / 4 ) (44) 

Equivalent results are obtained for the form factor considering only the contri- 
butions of the orbits involved in the bifurcation. 

In Fig. [14] we show the quantum results for Aoo versus energy E and for the 
form factor K(r) near ta, for the three values a = 9, 10 and 10.5 (as crosses, 
squares and circles, respectively). The solid line gives the universal GOE predic- 
tion, i.e., the first term in (14"0"1) . It agrees well with the quantum result at a = 9, 
in line with the near chaoticity of the system below the bifurcation. The dashed 
and dotted lines show the prediction (|44[) . which includes the bifurcating orbits 
A and L in the uniform approximation, and coincide well with the quantum 
results at and above the bifurcation. At the bifurcation (a = 0:4 = 10) where 
e = 0, Eq. (|44[) is consistent with the diagonal approximation for the bifurcating 
orbits and thus the same as that used in Ref. [13] • 




2000 4000_ 6000 8000 0.012 0.014 0.016 0.018 0.02 

E T 



Figure 14: Left: Saturation as a function of the energy. Right: Form factor near 
ta- Crosses, squares and circles mark quantum-mechanical results for a = 9, 10 
and 10.5, respectively. Bold line: GOE result, dashed and dotted lines: Eq. 



Figure[T4lmoreover shows that slightly above the bifurcation, i.e. at a = 10.5, 
the additional terms in Eq. (144[) . playing a role for e / 0, give a noticeable 
contribution, as seen by the dashed line. The main contribution comes from the 
term linear in e which is the nondiagonal contribution of the pairs of separate 
orbits A and L above the bifurcation. To see this, we evaluate their nondiagonal 
contribution in the Gutzwiller approximation for isolated orbits, which would 
become 

AL(„o„-dia g) _ 2 ^A/ Sin (*S\ V y/2 AS e 

°°( Gutz ) - TriTi \ Sm U ) ' ~ 2^\a\h ( ' 

(Although the diagonal contribution diverges at the bifurcation, the non-diagonal 
contribution stays finite there.) The additional factor l/v2> compared to the 
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last term in Eq. (|44[) . is due to the fact that the Gutzwiller approximation is 
not yet valid in this vicinity of the bifurcation (in particular, the difference in 
Maslov indices is different from the value 1 reached only far from the bifurcation 
where AS ^> % and hence e 3> 1). 

We see therefore that the nondiagonal contribution of the bifurcating orbits 
to the saturation value Aoo is non-negligible in a neigborhood above the bifur- 
cation. Note that the value of A^ is slightly enhanced also by the fact that the 
particular combination of Bessel functions in the uniform approximation (|36[) 
can be expressed by an Airy funtion (and its derivative, cf. Ref. [25]), which 
has its maximum slightly above the bifurcation. This effect is, however, not 
sufficient to explain the enhancement of A m found in our results, so that we 
can argue that the nondiagonal contribution is substantial. 

It is important to mention that this nondiagonal contribution exists as long 
as T% remains finite. In the strict semiclassical limit Ti — ► 0, the global uniform 
approximation (|36|) merges into the Gutzwiller trace formula for non-zero AS, 
and sm(AS/h) oscillates very fast, so that after the coarse-graining, the non- 
diagonal contribution will tend to zero. This is expected, since in the semiclas- 
sical approximation for mixed systems (Eq. [24]) , periodic orbits with different 
stability give rise to independent statistics. 

5 Conclusions 

In this case study we worked out for the quartic oscillator how (pitchfork) bi- 
furcations affect the density of states and thereby further measures of spectral 
correlations. This requires, at a first stage, detailed knowledge about the clas- 
sical bifurcation scenario in that system. At a second stage, we performed a 
comprehensive semiclassical calculation for the density of states invoking uni- 
form approximations for the bifurcating orbits involved. All features of the 
coarse-grained quantum density of states are adequately, and to high precision 
(mean level spacing), semiclassically reproduced, which is not evident in such 
a system with mixed phase space dynamics. Our semiclassical evaluation of 
the spectral rigidity close to the bifurcation shows strong deviations from the 
RMT behaviour, even though the phase space is predominantly chaotic and 
the bifurcation-affected phase space region appears negligible. This confirms 
that spectral statistics is rather susceptible with respect to bifurcation effects. 
Moreover we could unreval the role of orbit pairs born at the bifurcation which 
prevail with near-degenerate actions for larger control parameter regimes and 
strongly affect the spectral statistics. Such orbit pairs are obviously classically 
correlated and require a treatment beyond the diagonal approximation. 

This analysis moreover implies that in a comprehensive semiclassical ap- 
proach to spectral correlations in mixed systems, which still remains as a chal- 
lenge, off-diagonal contributions in the occuring multiple sums over periodic 
orbits should be considered, analogously to the purely hyperbolic case. 

Further open questions not answered in the present work include a corre- 
sponding analysis of how eigenstates are affected at a bifurcation. Finally, stud- 
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ies of bifurcation signatures in other observables such as quantum transport are 
still rare [15j and remain to be explored. 
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